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Novel materials and devices in photonics have the potential to revolutionize optical information 
processing, beyond conventional binary-logic approaches. Laser systems offer a rich repertoire of use¬ 
ful dynamical behaviors, including the excitable dynamics also found in the time-resolved “spiking” 
of neurons. Spiking reconciles the expressiveness and efficiency of analog processing with the robust¬ 
ness and scalability of digital processing. We demonstrate that graphene-coupled laser systems offer 
a unified low-level spike optical processing paradigm that goes well beyond previously studied laser 
dynamics. We show that this platform can simultaneously exhibit logic-level restoration, cascad- 
ability and input-output isolation—fundamental challenges in optical information processing. We 
also implement low-level spike-processing tasks that are critical for higher level processing: temporal 
pattern detection and stable recurrent memory. We study these properties in the context of a fiber 
laser system, but the addition of graphene leads to a number of advantages which stem from its 
unique properties, including high absorption and fast carrier relaxation. These could lead to sig¬ 
nificant speed and efficiency improvements in unconventional laser processing devices, and ongoing 
research on graphene microfabrication promises compatibility with integrated laser platforms. 


Recently, there has been a pertinacious exploration of 
the unifying boundaries between information communi¬ 
cation (dominated by optics) and information process¬ 
ing (dominated by electronics) in the same medium. In 
the context of information processing, nonlinear dynam¬ 
ical system^ 3 have been receiving considerable atten¬ 
tion due to their isomorphism to biological networks. 
Compared to binary-logic based methods implemented 
on standard von Neumann architectures, unconventional 
processing paradigms that are neuroinspired 3 ® are rela¬ 
tively more effective for solving certain tasks, such as pat¬ 
tern analysis, decision-making, optimization, and learn¬ 
ing. A sparse coding scheme, called spikingf 7 ® has re¬ 
cently been recognized by the neuroscience community 
as an important neural coding strategy for information 
processing. 8 ^ The continued evolution of photonic tech¬ 
nologies has reawakened interest for a relentless search 
in neuro-inspired optical information processing 2 12 1® to 
complement and enable new opportunities 16 17 ^ and po¬ 
tentially bridge the gap with information communication 
in the same substrate.® 

In this manuscript, we provide the first unified, exper¬ 
imental demonstration of low-level spike processing 7 8 
functions in an optical platform. We exploit uncon¬ 
ventional (excitable) dynamical properties of graphene 
laser systems to demonstrate the following features (in 
a fiber-based prot otype ) which are key impediments to 
optical computing® 18 : logic-level restoration, cascad- 
ability, and input-output isolation. Although a num¬ 
ber of approaches have demonstrated these properties 
separately,®*^ 1 no reported devices have simultaneously 
demonstrated these critical functionalities together in a 
single device!® 18 Our experimental prototype also pos¬ 
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sesses properties useful for processing tasks, including 
temporal integration, and sharp thresholding, leading to 
a very simple temporal classifier. 22 We include a simu¬ 
lation model that explains all of the observed behaviors: 
integration, thresholding, refractoriness, and pulse gen¬ 
eration. We also propose and simulate an analogous inte¬ 
grated device structure that exhibits the same dynamics 
in < mm 2 footprints. Scaling down the cavity length and 
overall size (by factors of millions) allows an integrated 
graphene excitable laser to exhibit dynamics on the or¬ 
der of ps timescales. Our model draws inspiration from 
novel insights in event-based information representation, 
dynamical excitability, and the unique material proper¬ 
ties of graphene. 

Spiking is a sparse coding scheme with firm code- 
theoretic justifications. 23 2 ^ Information is encoded in 
the temporal and spatial relationships between short 
pulses (or ‘spikes’). Spike codes—which are digital in 
amplitude but analog in time—exhibit the expressive¬ 
ness and efficiency of analog processing with the ro¬ 
bustness of digital communication. Spikes are typi¬ 
cally received and generated by nonlinear dynamical sys¬ 
tems, and can be represented and processed dynami¬ 
cally through excitability—a far-from-equilibrium non¬ 
linear dynamical mechanism underlying all-or-none re¬ 
sponses to small perturbations. 26 Excitable systems pos¬ 
sess unique regenerative properties and have been em¬ 
ployed for sensing microparticles with an optical torque 
wrench 27 and image processing utilizing a photosensi¬ 
tive Belousov-Zhabotinsky reaction. 28 In the context of 
spike processing, excitable laser system^® ® ® have 
been studied with the tools of bifurcation theory.® I® I®. 
Many dynamical systems that are explored are closely 
tied to underlying device physics, and, as such, the search 
for useful systems of this kind often involves novel mate¬ 
rials. 

Our approach exploits the unique properties of 
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FIG. 1. Graphene excitable fiber laser. The cavity consists of a chemically synthesized (see Methods) graphene SA 
(GSA) sandwiched between two fiber connectors with a fiber adapter and a 75-cm long highly doped erbium-doped fiber (EDF) 
as the gain medium. The EDF is pumped with a 980 nm laser diode via a 980/1550 nm wavelength-division multiplexer 
(WDM). An isolator (ISO) ensures unidirectional propagation. A polarization controller (PC) maintains a given polarization 
state, improving output pulse stability. The 20% port of an optical coupler provides the laser output at 1560 nm. To 
induce perturbations to the gain, 1480 nm excitatory pulses are incident on the system via a 1480/1550 nm WDM. These 
analog inputs—from other excitable lasers, for example—are directly modulated with an arbitrary waveform generator. Right: 
illustration of different possible phase-space dynamics associated with the laser intensity / and gain G of the system (see Text), 
as the various physical parameters (pump power, length of cavity, absorption) are varied. The desired excitable behavior, 
corresponding to the second phase-space schematic (red), is achieved when the parameter regimes drive the system toward a 
so-called homoclinic bifurcation. 1 


graphene, whose remarkable electrical and optical prop¬ 
erties have enabled several disruptive technologies. 32 34 
Graphene transistors are poised to be smaller and faster 
compared to their silicon counterparts P^ 36 but poor 
on/off current ratios resulting from a zero bandgap poses 
a serious challenge for conventional digital logic. Instead 
of using graphene’s electrical properties as an active ele¬ 
ment in conventional processing applications, we exploit 
its passive and unique optical properties to enable uncon¬ 
ventional processing. Since its emergence as a new type 
of saturable absorber (SA), graphene has been rigorously 
studied in the context of passive mode locking and Q- 
switchingj- 37 40 and has been preferred over the widely 
used semiconductor saturable absorbers 41 due to its high 
saturable absorption to volume ratio! 32 Graphene pos¬ 
sesses a number of other important advantages that are 
particularly useful in the context of processing, including 
a very fast response time, wideband frequency tunabil- 
ity (useful for wavelength division multiplexed networks), 
and a tunable modulation depth. Furthermore, graphene 
also has a high thermal conductivity and damage thresh¬ 
old compared to semiconductor absorbers. 

This work experimentally validates the theoreti¬ 
cally discovered^ SU dynamical isomorphism between 
semiconductor photocarriers and neuron biophysics, 
along with recent predictions of spike processing en¬ 
abled by graphene. 44 45 Ongoing research on graphene 
microfabrication 32 '"HI ESI ma y make it a standard tech¬ 
nology accessible in integrated laser platforms, which, to¬ 
gether with a suitable networking platform,L 4 2 could lead 


to a scalable platform for optical computingf^^ 


RESULTS 


Dynamical model. The dynamical system underly¬ 
ing the behavior of our spike processing unit is a gain- 
absorber cavity model, describing single mode lasers with 
gain and SA sections. Despite its simplicity, it can ex¬ 
hibit a large range of dynamical behaviors pO an d has 
been investigated in various contexts as the basis for 
an optical processor.^ The system, in its simplest form, 
can be described using the following undimensionalized 
equationP^ ® 


G(t) — qdd — G(t) — G(t)I(t)\ + 0(t) (la) 

Q(t) = lQ[B - Q(t) - aQ(t)I(t)\ (lb) 

Hi) = 7 i[G(t) - Q(t) - 1 ]I(t) + ef(G) (lc) 


where G(t) models the gain, Q(t) the absorption, and 
I(t) the laser intensity. A is the gain bias current, B is 
the absorption level, q <2 is the gain relaxation rate, qQ is 
the absorber relaxation rate, qj is the inverse photon life¬ 
time, and a is a differential absorption relative to the gain 
factor. We represent the spontaneous noise contribution 
to intensity via e/(G), for small e, and time-dependent 
input perturbations as 6(t). 

When the dynamics of pulse generation are fast com¬ 
pared to the dynamics of the gain medium, one can com¬ 
press the internal dynamics and obtain an instantaneous 
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FIG. 2. Proposed integrated graphene excitable laser, (a) Architecture of the hybrid InGaAsP-graphene-silicon 
evanescent laser (not to scale) with a terraced cutaway of the center. The device comprises a III-V epitaxial structure with 
multiple quantum well (MQW) region bonded to a low-loss silicon rib waveguide that rests on a silicon-on-insulator (SOI) 
substrate with a sandwiched heterostructure of two monolayer graphene sheets and an hexagonal boron nitride (hBN) spacer. 
Note that, unlike the fiber laser of Fig. ^ the gain section of this structure is electrically pumped. The cavity and waveguide 
are formed by the presence of a half-wavelength grating in the silicon, (b) Cross-sectional profile of the excitable laser with an 
overlaid electric field (E-field) intensity |E | 2 profile. The optical mode of the laser lies predominately in the silicon waveguide 
with a small portion of the mode overlapping the QWs of the III-V structure for optical gain and the 2D materials heterostructure 
for absorption. The silicon waveguide has a width, height, and rib etch depth of 1.5/xm, 500nm, and 300nm, respectively. The 
calculated overlap of the optical mode with the silicon waveguides is 0.558 while there is a 0.00046 overlap in the graphene 
sheets and 0.043 overlap in the QWs. The E-field intensities are calculated at a wavelength of 1.5/im. Graphene’s thickness of 
0.335nm and absorption coefficient of 301,655 cm -1 are used for the simulations (see Methods). 


pulse-generation modeP^l 

= G(G(t)-A) + 9(ty, (2a) 

if G(t) > Gthresh then (2b) 

release a pulse, and set G(t) G rese t 

where the input 6(t) can include spike inputs of the form 
6{t) = JA Si(t — Ti) for spike firing times t*, Gthresh is the 
gain threshold, and G reset ^ 0 is the gain at transparency. 

This system is analogous to a leaky integrate-and- 
fire (LIF) neuron model, commonly employed in com¬ 
putational neuroscience for modeling biological neural 
networks. Although it is one of the simpler spike- 
based models, the LIF model is capable of univer¬ 
sal computationsand the transmission of informa¬ 
tion through spike timings.^ 1 The gain-absorber system 
has been predicted to exhibit cascadability, logic-level 
restoration, and input-output isolation,^ satisfying the 
basic criteria for optical computing.^ 

Excitable Laser Systems. Our demonstration of spike 


processing is based on a graphene fiber ring laser plat¬ 
form (Fig. [I]). For comparison, we also perform numerical 
simulations of an analogous proposed integrated device 
(Fig. §■ Both devices along with their respective sim¬ 
ulation models and parameters are described in detail 
in the Methods section. The fiber ring laser contains an 
erbium doped fiber amplifier (gain section) and liquid ex¬ 
foliated graphene (absorber section), interacting with one 
another in a fiber ring (cavity). The ring laser pulses pe¬ 
riodically if driven above a threshold, modulated by the 
passive saturation of graphene absorption. This behavior 
has been studied in the context of high power, wideband 
passively Q-switched lasers, for which graphene has many 
favorable properties.^ 

The integrated device contains electrically pumped 
quantum wells (gain section), two sheets of graphene (ab¬ 
sorber section), and a distributed feedback-grating (sec¬ 
tion). In this design, we consider a hybrid silicon III-V 
laser platform in which the graphene layers are sand¬ 
wiched in between the silicon and III-V layers. The hy¬ 
brid III-V platform is highly scalable and amenable to 
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FIG. 3. Excitable dynamics of the graphene fiber laser. Note the blue and red curves correspond to input and output 
pulses, respectively, (a-c) Excitatory activity (temporal integration of nearby pulses) can push the gain above the threshold 
(measured to be ~ 275 nJ), releasing spikes. Depending on the input signal, the system can lead to (a) repeatable pulse 
generation, or suppressed response due to the presence of either (b) sub-threshold input energies (integrated power f \0(t)\ 2 dt) 
or (c) a refractory period during which the gain recovers to its resting value and the laser is unable to pulse (regardless of 
excitation strength), (d) Typical bursting behavior i.e. emission of doublets (two spikes) and triplets (three spikes) when a 
strong input drives the system over the threshold to fire repetitively, (e-f) Restorative properties (repeatable pulse reshaping) 
of spike processing; inputs with either the (e) same or (f) different energies. Test conditions: (a) three 4 /as input pulses 
separated by 3 /xs followed by a 10 /xs pulse after 52 /xs delay, and two 5.5 /xs pulses separated by 4.5 /xs after 50 /xs delay; (b) 
three 4 /xs input pulses separated by 3 /xs followed by a 4.5 /xs pulse after 52 /xs delay, and two 5.5 /xs pulses separated by 4.5 /xs 
after 50 /xs delay; (c) five 10 /xs input pulses separated by 10, 50, 20, and 25 /xs delays, respectively; (d) three input pulses 
with widths 25, 40, and 80 /xs separated by 55 /xs; in (a-d) the 980 nm pump is biased at 60.4 mW. (e) three input pulses with 
energies 252, 264, and 256 nJ separated by 85 and 65 /xs, result in outputs with pulses widths (i) 7.4, (ii) 6.1, and (iii) 6.9 /xs, 
respectively; 980 nm pump is biased at 61 mW. (f) three input pulses with energies of 265 nJ separated by 50 and 65 /xs, result 
in outputs with pulse widths (i) 6.3, (ii) 6, and (iii) 6.5 /xs; 980 nm pump is biased at 67.6 mA. 


both passive and active photonic integratiorP^. The in¬ 
tegrated device is capable of exhibiting the same behav¬ 
iors as the fiber prototype, but on a much faster time 
scale and with lower pulse energies. Figure ^l] compares 
the pulse repetition rate and pulse widths as a function 
of input power between the integrated device and fiber 
laser. In both cases, the rate of output pulses depends 
monotonically on the amount of power being consumed. 
This has many behavioral similarities with the behavior 
of rate neurons, which code information through spike 
frequency modulation.^ Although both lasers consume 
similar amounts of power, the integrated device pulses 
~ 10 6 times faster. This corresponds to a ~ 10 6 times 
decrease in the energy consumed per pulse. The devices 
(and their respective simulation models) are described in 
more detail in the Methods section. 

Excitability. We demonstrate that both the fiber ring 
laser and the integrated device are excitable and capable 
of performing spike processing tasks. Excitability is de¬ 


fined by three main criteria: (i) an unperturbed system 
rests at a stable equilibrium; (ii) a perturbation above 
the excitability threshold triggers a large excursion from 
this equilibrium; (iii) the system then settles back to the 
attractor in what is called the refractory period, after 
which the system can be excited again. 55 

Figures [3^i-c demonstrates excitability within the fiber 
ring laser. In this system, an excitatory pulse increases 
the carrier concentration within the gain region by an 
amount proportional to its energy (integrated power) 
through gain enhancement. Beyond some threshold ex¬ 
citation energy, the absorber is saturated, resulting in 
the release of a pulse. This is followed by a relative re¬ 
fractory period during which the arrival of a second ex¬ 
citatory pulse is unable to cause the laser to fire as the 
gain recovers. The system is also capable of emitting 
spike doublets or triplets (see Fig. [3 i) in which the inter¬ 
spike timing encodes information about the pulse width 
and amplitude, a useful encoding scheme for selective 
activation! 56 
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FIG. 4. Second-order properties of excitability, (a) Response of the excitable laser (output pulse energy) to a second 
input pulse as a function of the interspike interval between two identical (first and second) excitatory pulses after the first pulse 
has triggered an excitable response. Both the fiber based excitable laser experimental results and the integrated excitable laser 
simulations exhibit absolute refractory periods (where the second pulse produces no output) and relative refractory periods 
(where the output response is reduced from its resting value). However, the later operates ~ 10 3 faster (ns compared to fis) with 
~ 10 3 lower output pulse energies (pJ compared to nJ). Insets show the transient dynamics of the integrated excitable laser, i.e. 
the intensity I(t) and recovery of gain carriers G{t) as a result of input signal 0(t ), before (i) and after (ii) its refractory period, 
(b) Excitable laser’s output behavior in response to a single input pulse with different energies. The integrated excitable laser 
simulations also follows a similar relationship profile but with output pulse widths ~ 10 6 smaller (ps compared to fi s). (c) 
Response of the excitable laser when implemented a coincidence detector: the excitable laser is biased such that it will not fire 
unless two excitatory pulses are temporally close together. Output response is strongly dependent on the temporal correlation 
of two inputs. Average input power is kept constant with changing pulse interval. Insets show the simulated pulse dynamics 
for the integrated laser for pulses that are (i) closer vs (ii) farther apart. 


Since excitable systems are self-triggered, they exhibit 
important restorative properties. Different input pertur¬ 
bations often result in the same output, an important 
criteria for cascadability. Figures [3^ and f illustrate the 
response of the device as a result of a variety of input 
pulses. The excitable system responds in a stereotyped 
and repeatable way; all emitted pulses having identical 


pulse profiles. Outputs trigger asynchronously from in¬ 
put pulses, preserving analog timing information. 

Figure [4] show some key behaviors associated with ex¬ 
citability. Fig. [4^i provides information about the refrac¬ 
tory period for both the fiber laser and integrated device, 
which sets an upper bound on the pulse rate for a given 
unit. Similarly, Fig. [4]3 shows the output pulse width as 
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FIG. 5. Temporal pattern recognition, (a) Simple circuit with two cascaded graphene excitable lasers, (b) Measured 
output pulse peak power and pulse duration as a function of the time interval between the two input pulses, (c) Measured 
input and output waveforms at specific instances: (i) A t — r — —45/xs, (ii) A t ^ r — 135 fis, and (iii) At — r — 35 ns. The 
output pulse energy is largest when At « r showing the system only reacts to a specific spatiotemporal input pattern. 


a function of an input pulse for both the integrated and 
fiber lasers. The fiber experiment is corroborated with 
matching simulation results (see Methods). Although the 
pulse profile stays the same, its amplitude may change de¬ 
pending on the value of the perturbation. The integrated 
device exhibits the same behavior on the much faster time 
scale, recovering in ^ 1 nanosecond with pulse widths in 
/is, a factor of about ^ 10 6 faster than the fiber prototype 
in both respects. 

Temporal pulse correlation is an important processing 
function that emerges from excitability. A integrating ex¬ 
citable system is able to sum together multiple inputs if 
they are close enough to one another in time. This allows 
for the detection of pulse clusters, or potentially, coinci¬ 
dence detection of pulses across channels through the use 
of incoherent optical summing. 4 ^ Coincidence detection 
underlies a number of processing tasks, including asso¬ 
ciative memoryp 7 and a form of temporal learning called 
spike timing-dependent plasticity (STDP)P^^ Tempo¬ 
ral pulse correlation in the fiber laser experiment and 
simulation, and integrated laser simulation are shown in 
Fig. |4jj. Reducing the time interval between input pulses 
(i.e. simultaneous arrival) results in an output pulse. Al¬ 
though the fiber laser can function at kHz speeds, the 
internal dynamics of the integrated device allow it to 


function much faster, putting it in the GHz regime. 

Temporal pattern recognition. We demonstrate a 
simple pattern recognition circuit using several intercon¬ 
nected graphene fiber lasers. Pattern recognition of spa¬ 
tiotemporal phenomena is critical in the real-time pro¬ 
cessing of analog data. In the context of biological neu¬ 
ral systems, networks of spiking neurons convert ana¬ 
log data into a spikes and recognize spatiotemporal bit 
patterns.^ Spatiotemporal patterns play an important 
role in both visual 61 and audio 62 functionalities, and un- 
derly the formation of polychronized groups in the con¬ 
text of learning. 63 

As shown in Fig. IT’ we construct a simple two-unit 
pattern recognition circuit by cascading two excitable 
graphene (sic.) lasers with a delay r between them. In 
our case the objective is to distinguish (i.e. recognize) 
a specific input pattern: a pair of pulses separated by 
a time interval At ~ r, equal to the delay between the 
excitable lasers. Coincidence detection provides the dis¬ 
criminatory power for classification. 

A pulse doublet travels to both lasers, created using a 
modulator and arbitrary waveform generator. The out¬ 
put from the first laser travels to the second through a 
long single mode fiber (^km), which acts as a delay ele- 
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ment. The second laser is biased with a larger threshold 
such that it will not fire unless two excitatory pulses— 
the original input and output from the first laser—arrive 
at the same time (At ~ r = 135/is) (see Fig. [5J)). Syn¬ 
chronous arrival of these two spikes causes the release of 
a pulse. Experimental time traces for the inputs and out¬ 
puts are shown in Fig. A well-formed output pulse 
appears only for the desired two-pulse pattern. 

One can reduce the occurrence of pulses that are at 
a non-normalized amplitude through a sharper thresh¬ 
old function (i.e. Fig. [5 Jd) . This ratio can be optimized 
for application-specific purposes. The first laser acts as a 
nonlinear stage to simply regenerate input pulses because 
it is biased close to the excitable threshold. The second 
laser, on the other hand, requires two coincident pulses 
to reach its threshold. It therefore plays the role of a pat¬ 
tern classifier. Between the laser stages, a photodetector 
(PD), rather than direct optical input, modulates the 
laser driver (allowing wavelength conversion from 1560 
to 1480 nm). This PD-driven architecture (see outlined 
dashed box in Fig. has been explored in an inte¬ 
grated contexfP^, as a potential route to scalable on-chip 
networking. 47 ^ The dynamics introduced by the PD are 
analogous to synaptic dynamics governing the concentra¬ 
tion of neurotransmitters in between signaling biological 
neurons.^ This simple circuit demonstrates several im¬ 
portant features necessary for robust optical processing: 
well isolated input/output ports allow for the construc¬ 
tion of feedforward networks, and the spatio-temporal 
recognition of spikes allows the system to classify pat¬ 
terns. More complex recognition and decoding would be 
possible as the system is scaled. 

Stable recurrent circuit. We also demonstrate a self¬ 
recurrent graphene laser that can sustain a pulse trav¬ 
eling around the loop ad infinitum, providing a proof- 
of-principle demonstration of cascadability and pulse re¬ 
generation. Recurrently-connected dynamical networks 
which evolve toward a stable pattern over time (i.e. at¬ 
tractor networks) can exhibit hysteresis, and play a crit¬ 
ical role in memory formation and recall.^51 Equivalently, 
since a single unit with a self-referent connection can be 
mapped to an infinite chain of lasers, this system can 
be viewed as a demonstration of stability in arbitrarily 
many layered feed-forward networks. 

Figure [6^t illustrates an excitable graphene laser with 
a self-referent connection. The output is fed back to the 
input via single-mode fiber which acts as a delay element 
(100 /is) with an electronic weight W controls the modu¬ 
lation depth of the PD providing an all-or-none response 
depending on whether it is above or below a given thresh¬ 
old. Figure depicts the system’s ability to demon¬ 
strate bistability when feedback is present. It is capable 
of settling to an attractor in which a single pulse trav¬ 
els around the loop indefinitely. This circuit represents a 
test of the network’s ability to handle recursive feedback, 
and the stability of the pulse is a sign that the system is 
cascadable. 



Time (|js) 


FIG. 6. Self-recurrent bistable circuit with the 
graphene excitable laser, (a) Setup to test the self¬ 
referent connection, (b) Input and output waveforms. The 
first output pulse is fed back to the input after being delayed 
by ~100 /is, which initiates another excitatory pulse at the 
output. This recursive process results in a train of output 
pulses ad eternum at fixed intervals. Inset shows an output 
pulse profile and sech 2 fitting curve. 


DISCUSSION 


We have demonstrated that the complex dynamics of 
graphene excitable lasers can form a fundamental build¬ 
ing block for spike information processing. In addition 
to single-laser excitability, we showed two instances of 
key spike processing circuits: temporal pattern recogni¬ 
tion and stable recurrence. A photonic coincidence de¬ 
tection circuit forms the building block of the spatiotem- 
poral pattern recognition circuit we have also demon¬ 
strated by cascading two excitable lasers as computa¬ 
tional primitives. This simple demonstration of tempo¬ 
ral logic implies that spiking neural networks of such 
excitable lasers are capable of categorization and deci¬ 
sion making. Combined with learning algorithms such as 
STDP, networks could potentially perform more complex 
tasks such as spike-pattern cluster analysis.^ 3 A bistable 
recurrent spiking circuit enabled by the graphene ex¬ 
citable laser shows that processing networks of excitable 
lasers are capable of indefinite cascadability and informa¬ 
tion retention, a pre-requisite for more complex types of 
temporal attractors in recurrent networks. In networks 
of more lasers, spiking attractors can be more numer¬ 
ous, complex, and even competitive in order to achieve 
different information processing goals. 

Ongoing research on graphene microfabrication could 
make it a standard technology accessible in integrated 
platforms. We proposed an integrated graphene- 
embedded cavity design and adapted the fiber model of 
excitability to a semiconductor device model. Our re¬ 
sults show that an integrated device could maintain the 
essential behaviors required for spike information pro¬ 
cessing while reaping significant energy and speed im- 








































provements, potentially opening up applications for bio¬ 
logically inspired adaptive algorithms in presently inac¬ 
cessible regimes of computings 


METHODS 

Fiber laser simulation. To simulate the fiber laser, we 
constructed rate equations based on the carrier dynamics 
in an EDF amplifier, roundtrip intensity, and loss. The 
dynamics of an EDF can be described using the following 
equations for fractional excited state population 77 , 2 , frac¬ 
tional ground state population m, and k optical beams 
of intensities 7/J®!: 


dn 2 

dt 
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&aklk _ 
HuJi 


-ni 


(M) - E 


&eklk 

hu k 


n 2 {z,t) 


n 2 jt) 

r n 


(3) 


Each term represents the transition rate for each photon, 
where where h is Planck’s constant and ujk is the fre¬ 
quency of mode k. cr a k and a e k represent the absorption 
and emission cross sections of each mode fc, respectively, 
and the fractional populations satisfy ni = 1 — n 2 . Our 
interest is in the modes at pump wavelengths 980 nm and 
1480 and lasing modes which hover around 1520-1530 
nm. We define pump intensity I p at 980 nm, and input 
signal intensity I s at 1480 nm, and the round trip inten¬ 
sity I r at 1550 nm. Although the fiber laser is largely 
multi-mode, the modes are closely spaced to one another 
and possess similar cross sections. We can therefore ap¬ 
proximate these modes with a single roundtrip intensity 
1= Ik equal to the sum of lasing modes, and define ef¬ 
fective cross sections a er and a ar . We also use a lumped 
approximation and represent the carrier density as a sin¬ 
gle variable!®" We average over the fiber length z to arrive 
at the following differential equation for average carrier 
density = f n(z)dz: 


dn 2 

dt 


+ Wa S ni(t) - a es n 2 {t)]T- 
riujp nw s 

r _ , s _ , s-| I r 77*2 (t) 

+ Warni(t) - cr er n 2 (t)\ — - 

llLUry. T n 


(4) 


Powers Pk injected into the erbium-doped section are 
related to the average intensity within this section via 
4 = %[(e Sfc(t) - 1 )/A e f f g k (t)]Pk(t), where g k (t) « 
r/ c rit[cr e / c n 2 (t) — (T a kIpif)\LEr is the gain experienced by 
the mode over the length of the erbium fiber, A e ff is 
the effective cross sectional area of the fiber, rjk is the in¬ 
jection efficiency, T^ is the confinement factor, n t is the 
erbium ion density, and Le t is the length of the erbium 
section. 67 

We can define roundtrip equations for both round-trip 
loss q(t) and round-trip power averaged over the fiber 


TABLE I. Fiber Laser Param. (from data and [38] 1 


Param. 

Description 

Value 

Ip 

pump intensity 

4.39 x 10 8 W/m 2 

Pk 

pump coupling efficiency 

.21 

Vs, V p 

input sig. and pump freq. 

194, 306 THz 

T n 

erbium lifetime 

9 ms 

C ap 

abs. cross section (pump) 

2.87 x 10“ 25 m 2 

O'asi <Tes 

abs., em. cross sec. (sig.) 

3.01, .948 xlO -25 m 2 

Tar 1 T er 

abs., em. cross sec. (rt.) 

3.21, 4.54 xlO -25 m 2 

^eff(r) 

eff. cross sec. area (rt.) 

2.77 x 10 -11 m 2 

Efe(p) 5 rfe(s) 

conf. factors (pump, sig.) 

.849, .638 

n t 

erbium ion density 

5.8 x 10 25 m 3 

LEr 

erbium fiber length 

75 cm 

qo 

small-sig. SA absorption 

.5 

T q 

SA lifetime 

2 ps 

Psat 

SA saturation energy 

10 pJ 

Tr 

cavity roundtrip time 

90 ns 

l 

cavity intrinsic loss 

1.1 

Psp 

spont. noise term 

5 W/s 


length P r (t): 

dq = q(t) - qp _ q(t)P r {t) 

dt Tq Psat 

dP 

Tr= [grit) - q(t) - l]P r (t) + Psp ( 6 ) 

where go represents the small-signal absorption of the 
SA, r q the absorber relaxation time, E sat the saturation 
energy, Tr the round-trip cavity time, l the round-trip 
loss, p sp a small spontaneous noise term, and g r (t) = 
r r n t [a er n 2 {t) — cr ar ni{t)]L e t the erbium fiber round-trip 
gain. Equations (|4]-[6| represented the model used for the 
simulation. Parameters are shown in Table |IJ These equa¬ 
tions were stepped iterately using Runge-Kutta methods 
to generate time traces and measure various properties. 

It is possible to recover the simplified, undimension- 
alized model that underlies the observed behaviors by 
noting that n 2 (t) does not change significantly over time 
(i.e. n 2 {t) = uk + 5n(t) where Sn(t) « uk) and sub¬ 
stituting g r (t) into equation ©• These approximations 
lead to a bilinear set of equations that are analogous to 
equations ( [Ta}jTc| ). 

Integrated-device simulation. Design principles are 
chosen for compatibility with recent graphene deposi¬ 
tion and patterning techniques.^ ^ The device epitaxial 
layer structure includes both quantum wells (QWs) and 
graphene coupled to a single optical mode, shown in Fig¬ 
ure [2] Both the QWs and graphene provide complemen¬ 
tary properties—whereas QWs provide high efficiency 
gain, the graphene provides strong, fast and wideband 
saturable absorption. The difficulty in coupling graphene 
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directly to the optical mode could be resolved for in¬ 
stance by wafer bonding a III-V laser on top of deposited 
graphene, avoiding any interaction between graphene and 
electrical pumping, as shown. For improved dynamics, 
we consider two pristine layers of graphene, protected by 
a atomically flat layers of boron nitride (BN) to prevent 
each graphene sheet from interacting too strongly with 
surrounding materials. We computed the optical mode 
of this structure using an eigenmode expansion (EME) 
technique. 

Using the confinement factor from above and other var¬ 
ious parameters, we simulated the device using a lumped 
rate equation model. Beginning with the theory for 
graphene, the behavior is well approximated by a sim¬ 
ple saturation model, given b}® 

Q^a) = S / - \~CtNS (7) 

where a(v a ) represents the absorption coefficient (per 
unit length), v a is the two-dimensional carrier density 
in graphene, v s the 2D saturation carrier density, a^s 
the saturable absorption, and a^s the non-saturable ab¬ 
sorption. The resulting rate equations are given by: 


= v g g(n g )N ph - v g a(v a )N ph + — + R sp ( 8 a) 
at T p h 


dn g _ I g + <j>{t) 


. Nph 71 


dt 


eV a 


--*+*(*) ( 8 b) 

V a T n 


dl^a Va / \ Nph 


( 8 c) 


N p h represents the number of photons in the cavity and 
rig the carrier density QW gain region. (Note: the 
variable v a represents the surface carrier density within 
graphene, chosen for convenience as graphene is two- 
dimensional.) g(n g ) and a(n a ) describe the gain and 
absorption per unit length, v g the group velocity, r the 
lifetimes, I g current pumped into the gain region, R sp a 
small spontaneous noise term, V g the volume of the gain 
region, A a the area of the graphene sheet, and (f>(t) an in¬ 
put current modulation term. The input power P g that 
is driving the laser can be computed by P g = I g x vl 
where vl is the voltage applied across the gain section of 
the laser. Gain and loss are assumed to take the forms: 


g(n g ) = T g g 0 log{n g /n tr ) 


a{v a ) 


0‘S 

1 + v<x!v s 


where T^ is the gain confinement factor, n tr is the trans¬ 
parency density (cm -3 ) and v s is graphene transparency 
density in two dimensions (cm -2 ). Non-saturable ab¬ 
sorption Oi^s is not included as it manifests as cavity 
losses, becoming absorbed into the photon lifetime r p h- 
Parameters are shown in Table [TT1 We simulated the rate 
equation model using Runge-Kutta methods. 


TABLE II. Hybrid Integrated Laser Param. (from simula¬ 
tions, and [53 , 691 1 


Param. 

Description 

Value 

A 

lasing wavelength 

1550 nm 

v g 

group velocity 

c/3.5 

Vg 

gain region volume 

2.55 xlO -12 cm 3 

A(x 

graphene sheet area 

1.5xl0 -6 cm 2 

r * 

gain region confinement factor 

0.034 

T g 

gain region carrier lifetime 

1.1 ns 

Tex 

graphene carrier lifetime 

405 fs 

T P 

photon lifetime 

2.4 ps 

go 

log gain coefficient 

972 cm -1 

Qt s 

waveguide saturable absorption 

150 cm -1 

Tltr 

3D gain region transparency density 1.75xl0 18 cm 3 

Vs 

2D graphene transparency density 

1.06xl0 13 cm” 2 

Rsp 

spontaneous noise term 

1 x 10 10 s” 1 

VL 

applied voltage (gain section) 

1.1 V 


One can recover the undimensionalized equations with 
several approximations and variable substitutions. Mak¬ 
ing a linear approximations to both the gain and ab¬ 
sorption, i.e. g{n g ) ~ go[n g — n t r\/n t r and a(u a ) ~ 
a£[l — Vu! z/ s ], leads to an equation analogous to the sim¬ 
plified gain-absorber model described by equations ([la]- 

0 - 

Excitable fiber ring laser cavity. The EDF employed 
in the laser cavity is a gain fiber (LIEKKI Er80-4/125), 
with peak core absorption coefficients of 60, 50 and 110 
dBm -1 at 980, 1480, and 1530 nm, respectively. It has a 
large area core with a mode field diameter of 6.5 gm at 
1550 nm and a core numerical aperture of 0.2. The high 
erbium ion doping concentration reduces the required 
fiber length significantly while providing strong gain and 
reducing nonlinear effects (four-wave mixing, stimulated 
Raman scattering, stimulated Brillouin scattering). The 
length of the EDF (75 cm) is chosen to ensure population 
inversion with the desired pump power so that the EDF 
does not play any role as a SA to realize excitability. All 
fibers used in the cavity are polarization-independent. A 
polarization controller consisting of three spools of SMF- 
28 fibers acting as retarders is used to maintain a given 
polarization state after each round trip improving the 
output pulse stability.^ 

Graphene sample preparation. Graphene samples 
are prepared by chemical reduction of graphene oxide 
(GO) with hydroxylamine hydrochloride (NH 2 OH 3 HCI) 
with a slightly modified recipe. 70 ! 25 ml of 0.5 mg/ml 
GO (Graphene supermarket #SKU-GO-W-175) is di¬ 
luted with 25 ml of deionized (DI) water, 200 gL of 28 
wt.% ammonium hydroxide (Sigma-Aldrich #338818), 
and 25 mg of NH 2 0H 3 HC1 (Sigma-Aldrich #431362) 
in a 100 ml round-bottom flask and stirred. The mix¬ 
ture is transferred to a water bath and heated at ^90° C 
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with stirring at 350 rpm for 90 mins. The color of the 
mixture changes from yellowish brown to homogeneous 
black and precipitating from solution indicating reduc¬ 
tion has taken place. The reduced GO (rGO) is filtered 
and washed three times with DI water. The rGO is then 
suspended in 50 ml of DI water with 50 mg of sodium de- 
oxycholate (Sigma-Aldrich 7^30970), and stirred until the 
salt dissolves. The mixture is sonicated for 30 mins using 
ultrasound-assisted functionalization resulting in a stable 
rGO suspension. A micropipette is used to transfer 5 fiL 
of the as prepared liquid to an angle-polished fiber con¬ 
nector (FC/APC). The deposited sample is dried with a 
heat gun operating at around 120°C for ^3 to 5 mins. 
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Figure S 1. Typical characteristics of the passively Q- 
switched fiber and integrated lasers, (a) Output pulse rep¬ 
etition rate and (b) pulse width as a function of the pump 
current. 








